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Abstract:  Adult  skeletal  muscle  is  the  dominant  system  through  which  complex  physical  actions 
are  accomplished  and  is  uniquely  capable  of  repair  and  regeneration  after  different  types  of  insult 
or  injury  via  resident  stem  cells.  The  in-vivo  epigenomic  and  transcriptional  mechanisms  through 
which  skeletal  muscle  repairs  itself  are  partly  understood  and  herein,  we  administer  severe 
muscle  trauma  to  a  mouse  tibialis  anterior  muscle  and  monitor  the  in-vivo  dynamics  of  three 
histone  modifications  and  gene  expression  across  nine  time  points  after  injury.  Integrating  the 
genome- wide  datasets,  we  show  extensive  chromatin  state  remodeling  at  c/.s-rcgulatory  elements 
and  observe  the  cooperative  activity  of  lineage-specific  transcription  factors  during  progressions 
through  different  stages  of  healing.  We  highlight  how  one  of  the  major  regenerative  programs 
(IGF  and  PI3K  signaling)  activated  after  injury  is  titrated  by  negative  regulators  temporally  on 
multiple  levels.  Fastly,  we  illustrate  the  dynamics  of  non-coding  RNAs  during  various  stages 
after  trauma,  most  notably  during  the  switch  from  myoblast  proliferation  towards  differentiation 
and  eventual  muscle  regeneration.  These  results  provide  an  unbiased,  comprehensive  view  of  the 
central  factors  that  regulate  muscle  regeneration  after  severe  trauma  and  underscore  the  multiple 
levels  through  which  both  local  and  global  genetic  and  epigenetic  patterns  are  regulated  to  enact 
appropriate  and  timely  repair  and  regeneration. 


Introduction 

Adult  skeletal  muscle  is  a  post-mitotic  organ  that  coordinates  movement  and  constantly  grows 
and  adapts  by  remodeling  its  structure  and  metabolism.  After  insult  or  injury,  adult  skeletal 
muscle  enables  repair  and  regeneration  of  existing  fibers  through  a  population  of  stem  cells  that 
reside  underneath  the  basal  lamina  called  satellite  cells1  (SCs).  The  SCs  reside  in  a  specialized 
niche2,2  and  change  their  quiescent  complexion  after  injury  via  various  environmental  signals  and 
communication  with  infiltrating  immune  cells4'5  and  fibro-adipocyte  progenitors  (FAPs)6.  Many 
essential  SC  regulators  have  been  studied  during  the  muscle  repair  and  regeneration  process, 
providing  an  excellent  model  to  study  the  actions  of  different  factors  that  support  myogenic 
repair  and  regeneration  sub-processes. 

It  is  well  established  that  multiple  repair  and  regeneration  sub-processes  accomplished  by  SCs 
(and  other  cell  types)  after  muscle  injury  are  orchestrated  by  distinct  transcriptional  networks7, 
encompassing  epigenetic89,  transcriptional10  and  post-transcriptional  events.  However,  the 
integrative  dynamics  of  transcriptional  networks  and  regulatory  epigenetic  switches  at  genome¬ 
wide  levels  have  not  been  characterized  in-vivo  and  as  such,  our  understanding  of  the  molecular 
processes  and  transcription  factors  (TFs)  involved  in  myogenesis  and  muscle  regeneration  have 
been  limited.  Profiling  the  expression  and  chromatin  state  of  ds-regulatory  elements11  after 
severe  muscle  trauma  provides  a  powerful  method  to  understand  the  molecular  determinants  of 
cellular  fate,  healing  progression  after  severe  muscle  trauma  and  can  provide  crucial  insights  for 
development  of  therapeutic  modalities  for  effective  repair  and  regeneration  of  severe  muscle 
trauma,  myopathies  such  as  Duchenne  muscular  dystrophy  and  aging. 

The  ds-regulatory  networks  that  orchestrate  in-vivo  muscle  repair  and  regeneration  after 
traumatic  injury  have  only  been  partly  characterized 1215  and  herein,  the  dynamic  in-vivo 
evolution  of  expression  and  three  different  chromatin  modifications  (H3K4me3,  H3K4mel  and 
H3K27ac)  were  profiled  across  nine  time  points  (t  =  3  hrs  to  672  hrs)  from  an  injured  and 
uninjured  contralateral  tibialis  anterior  muscle.  The  generated  genomic  maps  were  then 
contrasted  against  MyoD  and  MyoG  ChIP-seq  data  from  C2C12  cells14'15  to  determine  shared 
and  distinguishing  signatures  at  ds-regulatory  elements  during  different  stages  after  injury.  The 
dynamic  chromatin  state  transitions,  differential  binding  at  TF  motifs,  levels  of  numerous  coding 
and  noncoding  transcripts  and  corresponding  adaptive  molecular  responses  were  integrated  and 
assessed  to  construct  a  comprehensive  view  of  the  key  transcriptional  and  chromatin  factors  that 
influence  and  modulate  in-vivo  muscle  repair  and  regeneration  dynamics. 

Results 

Severe  muscle  trauma  induces  extensive  transcriptional  changes  and  chromatin  state  remodeling 

Figure  1  shows  chromatin  immunoprecipitation  followed  by  DNA  sequencing  (ChIP-Seq)  was 
used  to  globally  map  the  chromatin  state  of  various  c/.v-regulatory  elements  at  nine  time  points 
that  were  lumped  into  three  key  stages  after  a  traumatic  muscle  injury  (early:  3-24  hours  after 
injury,  middle:  48-168  hours  after  injury,  late:  336-672  hours  after  injury).  Validated  antibodies 
(Supp.  Fig.  1)  for  histone  H3  lysine  4  trimethylation  (H3K4me3),  a  modification  associated  with 
promoters,  H3K4mel  and  H3K27ac,  associated  with  poised  and  active  enhancer  regions16, 
respectively,  were  used  to  enrich  chromatin.  Each  tissue  from  each  time  point  was  imrnuno- 
precipitated  using  the  three  antibodies,  and  ChIP-Seq  maps  from  the  same  antibody  and  time 
point  were  merged,  resulting  in  54  chromatin-state  maps  covering  >1.5  billion  reads.  RNA-Seq 


2 


maps  for  both  the  injured  and  contralateral  tissues  for  each  time  point  were  also  generated  and 
analyzed17.  Hierarchical  clustering  of  the  RNA-Seq  data  through  time  and  peak  fold  change 
revealed  clusters  upregulated  at  different  time  periods  that  were  associated  with  different  stages 
of  muscle  repair  and  regeneration  (Fig.  lb).  Chemokine  ligands  CCL2  and  CCL7,  which  are 
important  for  the  recruitment  of  various  immune  cells  to  the  injured  muscle,  peaked  in 
expression  in  the  early  period  (at  24  h  after  injury).  Annexins  1  and  2  (Anxal  and  Anxa2), 
cellular  membrane  binding  proteins,  demonstrated  a  different  expression  profile  and  peaked  in 
the  middle  period  (between  48h  and  72h)  after  injury  suggesting  the  beginning  of  the 
regenerative  program.  Tropomyosins  3  and  4  (Tpm3  and  Tpm4),  thin  filament  proteins  that 
participate  in  muscle  contraction,  demonstrated  a  similar  temporal  expression  profile  to  the 
annexin  family  and  other  sarcomere  proteins,  which  peaked  in  the  middle  period  and  remained 
upregulated  into  the  late  stage,  in  line  with  a  productive  healing  process. 

Figure  2  shows  in  the  early  time  period,  93,149  sites  were  enriched  for  H3K4me3,  and  104,890 
sites  for  the  middle  period  followed  by  141,948  sites  for  the  late  period.  Of  these,  1,606  were 
differentially  enriched  in  the  early  time  points  relative  to  the  controls,  6,733  in  the  middle  time 
points,  and  6,121  in  the  late  time  points  (Supp.  Fig.  2).  To  integrate  sites  that  gained  or  lost  the 
H3K4me3  modification  during  the  time  course  with  transcriptional  activity  (FPKM>1),  RNA 
sequencing  (RNA-Seq)  results  were  performed  at  the  same  time  points17  and  contrasted  against 
the  genome- wide  maps.  Approximately  58%  percent  of  sites  were  found  to  associate  with 
transcriptional  activity  (FPKM>1)  and  4,241  sites  exhibited  dynamics  (acquisition  or  loss  of  the 
histone  modification)  through  at  least  one  stage  (early-middle,  middle-late,  early-late). 

In  contrast  to  H3K4me3,  which  has  previously  been  shown  to  be  largely  static  during  chromatin 
remodeling  events18,  H3K4mel  and  H3K27ac  demarcate  enhancer  elements,  which  are  highly 
dynamic19.  H3K4mel  resides  on  both  poised  and  active  enhancers,  as  to  where  H3K27ac  marks 
active  enhancers20.  In  the  early  time  period,  93,938  sites  enriched  for  H3K4mel  were  identified, 
106,353  sites  for  the  middle  period,  and  28,663  sites  for  the  late  period.  The  active  enhancer 
mark  H3K27ac  showed  similar  behavior  with  32,285  sites  enriched  in  the  early  time  period, 
43,780  sites  in  the  middle  period  and  13,748  sites  for  the  late  period.  Figure  2  depicts  the 
highest  number  of  sites  that  acquired  H3K27Ac  was  found  in  the  72hr  period  (28,178  sites) 
compared  to  1 1,535  for  the  48hr  stage  and  4,067  for  the  168hr  stage. 

Immune  Cell  Programs  Contribute  to  Differential  Chromatin  States  After  Severe  Trauma 

In  the  early  period,  numerous  transcripts  associated  with  inflammation,  invading  immune  cells4, 
cytokine  signaling,  and  apoptosis  were  detected.  Gene  Ontology  (GO)  tenn  analysis  of  the 
enriched  chromatin  peaks  in  the  early  period  that  underwent  changes  for  the  injured  samples 
relative  to  uninjured  controls  demonstrated  enrichments  similar  to  the  transcriptional  groups 
associated  with  immune  response,  chemokine  and  cytokine  signaling,  inflammation,  and 
metabolism  of  lipids  and  lipoproteins  (Fig.  3a  &  Supp.  Fig.  3).  The  increase  in  number  of 
enriched  chromatin  peaks  associated  with  immune  response  and  inflammatory  pathways  can 
partially  be  attributed  to  infiltration  of  monocytes,  lymphocytes  and  macrophages,  which  are 
essential  for  functional  recovery  and  generally  follows  a  two-stage  process.  In  the  first  stage,  the 
invading  immune  cells  induce  a  protective,  pro-inflammatory  state  and  after  several  days,  switch 
to  an  anti-inflammatory  program5.  In  line  with  this  observation,  transcription  factor  binding 
analysis  of  the  enriched  c/.s-regulatory  sites  showed  binding  motifs  for  cellular  stress  TFs  such  as 
NF-kP,  serum  response  factor  (SRF),  PU.l  and  AP-1  (Fig.  2d),  which  promote  chromatin 
accessibility  and  facilitate  additional  TFs  to  bind  and  modulate  expression  levels.  The  result  of 
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these  enrichments  at  c/.s-regulatory  elements  is  also  consistent  with  high  expression  of  IL-6  and 
TNFa,  both  of  which  induce  genes  such  as  c-Fos,  c-Jun,  Atfl,  Atf3  (Fig.  3c-d)  and  JunB  (Supp. 
Fig.  3).  Additional  binding  sites  were  also  observed  for  STATs  and  JAK-STAT  signaling,  which 
along  with  pro-inflammatory  cytokines  activate  downstream  target  genes  such  as  c-Myc, 
reported  to  drive  satellite  cell  proliferation21  and  repress  myogenic  TFs  such  as  MyoD  and  MyoG 
(Supp.  Fig.  4).  These  TF-coordinated  programs  most  likely  permit  the  injured  site  to  activate  and 
simultaneously  titrate  satellite  cell  expansion  and  immune  cell  infiltration  as  well  as  prevent 
premature  myogenic  differentiation. 

Tissue  Architecture  Cues  Influence  Chromatin  States  &  Transcription  Factor  Binding 

In  the  middle  time  period  (48-72hrs  after  injury),  upregulation  of  SMADs  and  TGF-(3  signaling 
was  observed  along  with  increases  in  expression  of  TFs  such  as  Tead4  and  MyoD.  These  two 
elements  act  to  repress  c-Myc  expression  and  promote  a  transition  from  proliferation  and  Notch 
signaling  towards  Hippo  and  Wnt  signaling22'23  (Fig.  4a).  SMADs  have  been  shown  to  interact 
with  chromatin  remodeling  complexes  such  as  histone  acetyltransferases  p300  and  CBP  (CREB- 
binding  protein)  to  induce  H3K27  acetylation,  which  is  consistent  with  the  observation  that  the 
highest  number  of  sites  that  acquired  H3K27ac  was  in  the  middle  time  period  (Fig.  Id). 

GO  analysis  of  the  differential  H3K27ac  peaks  in  the  injured  samples  for  the  middle  period 
revealed  multiple  enrichments  for  G-protein  coupled  receptors  (Fig.  4b)  and  the  Rho  family  of 
small  GTPases.  This  result  is  consistent  with  the  peak  in  expression  of  SRF  target  genes  and 
anti-inflammatory  cytokines  interleukin-4  (IL-4)  and  interleukin- 10  (IL-10),  which  induce  M2 
macrophage  polarization  and  are  essential  components  for  resolution  of  inflammation  and  tissue 
repair5  (Fig.  4a).  Intriguingly,  RhoA,  GTPase  and  SRF  each  selectively  modulate  the  expression 
of  MyoD24.  Increases  in  expression  of  these  factors  was  also  mirrored  by  upregulation  of  the 
macrophage-derived  matrix  metalloproteinase  12  (Mmpl2),  which  cleaves  and  inactivates  CXC 
chemokines  (Cxcll,  -2,  -3,  -5,  and  -8)  and  monocyte  chemotactic  proteins  (Ccl2,  -7,  -8,  and  -13), 
inhibiting  leukocyte  flux  to  the  injured  site  and  abrogating  the  amount  of  pro-inflammatory 
molecules  present  in  the  injured  tissue  (Supp.  Fig.  4).  Lastly,  differential  H3K27ac  enrichments 
(Fig.  4c)  were  viewed  for  genes  associated  with  extracellular  matrix  remodeling,  such  as 
components  of  the  basement  membrane  (myoferlin,  laminin,  collagen  VI  genes,  annexins),  and 
alter  the  stiffness  of  newly  formed  or  repaired  ECM.  Upregulation  of  GPCRs  and  altered  ECM 
stiffness25  have  previously  been  shown  to  stimulate  mechano-sensitive  Yap  and  Taz  to 
translocate  into  the  nucleus  and  induce  expression  of  genes  that  promote  satellite  cell 
proliferation26,  migration  and  Hippo  and  Wnt  signaling. 

The  collective  integration  of  these  different  enrichments  at  c/.s-rcgulatory  elements  for  the  middle 
period  suggests  the  elastic  nature  of  injured  microenvironment  (which  is  also  regulated  by 
invading  fibroblasts  and  different  types  of  immune  cells)  and  anti-inflammatory  soluble  signals 
provide  cues  for  nuclear  shuttling  of  TFs27  and  activation  of  enhancers  that  signal  the  resident  SC 
population  to  expand28  and  migrate  towards  damaged  fibers.  This  result  also  implies  that 
regenerating  muscle  possesses  a  plastic  chromatin  state  at  enhancers29,  which  can  provide 
flexibility  to  quickly  respond  to  different  ECM  cues  as  well  as  facilitate  sensitive  interactions 
with  invading  cell  types  that  may  further  influence  differential  usage  of  c/.s  -regulatory  elements. 

Transcription  Factor  Coordinated  Activation  of  Myogenic  Differentiation  and  Regeneration 

In  the  middle  period,  a  unique  and  complex  set  of  events  occurs  whereby  multiple  pathways 
begin  the  process  of  muscle  repair  and  regeneration.  Analysis  of  the  total  and  differential 
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H3K27ac  enrichments  for  the  middle  period  revealed  multiple  binding  sites  for  MyoD,  MyoG 
and  MEF2,  along  with  RunX,  Teads,  PPARy  and  various  ternary  complex  factors  (TCFs)  of  the 
ETS-domain  family.  At  enriched  promoter  peaks  for  both  the  MyoD  and  MyoG  datasets 
(K4me3  &  MyoD+  or  K4me3  &  MyoG+),  we  also  find  binding  motifs  for  Prdml4,  a  known 
histone  methyltransferase,  and  Pbx3,  a  TF  that  cooperatively  acts  with  MyoD  to  enact  fast- twitch 
muscle  expression30.  Next,  focusing  on  the  enhancer  sites  (H3K27ac+  &  H3K4me3‘)  during  the 
middle  period,  the  number  of  enriched  sites  that  overlapped  with  MyoD  binding  was  viewed  to 
peak  at  72hrs  (Fig.  5a).  This  observation  is  consistent  with  the  largest  increase  in  sites  that 
acquired  H3K27ac  was  at  72h  as  well  as  the  known  association  between  MyoD  and  p300/CBP 
on  E-box  motifs  of  target  genes  during  myogenic  differentiation.  The  increase  in  detected  MyoD 
binding  sites  at  enhancers  also  mirrored  the  temporal  activation  of  enriched  GO  terms  associated 
with  myogenic  differentiation  (Fig.  5c)  such  as  different  types  of  growth  factor  signaling  (insulin 
growth  factor  -  IGF,  fibroblast  growth  factor  -  FGF,  hedgehog),  p38  mitogen-activated  protein 
kinase  (p38-MAPK),  and  phosphatidylinositol  3-kinase  (PI3K)  signaling  (discussed  further 
below).  IGF  signaling  plays  a  crucial  role  in  muscle  repair  and  regeneration  after  injury  and 
previous  efforts  have  shown  that  activated  satellite  cells  express  IGF  binding  proteins  (IGFBPs) 
31.  Concordantly,  several  IGFBPs  and  IGF2  were  observed  to  be  upregulated  in  this  time  period 
(Fig.  5b).  The  p38  MAPK  pathway  has  also  previously  been  shown  to  stimulate  myogenic 
transcription  through  recruitment  of  the  chromatin-remodeling  complex  SWI-SNF  to  specific 
loci  and  neutralizes  TFs  such  as  c-Jun  that  inhibit  cell-cycle  exit32.  In  line  with  this,  we  view 
induction  of  cyclin-dependant  kinase  inhibitor  p2 1 ,  which  is  activated  by  MyoD  and  promotes 
cell-cycle  exit  and  GSEA  reactome  sets  associated  with  DNA  unwinding.  The  presence  of 
PPARy  enrichments  may  be  derived  from  resident  fibro/adipogenic  progenitors  that  expand  and 
proliferate  in  damaged  muscle  tissue  to  produce  factors  that  support  myogenic  differentiation6'33. 
The  various  pathways  converge  to  direct  TFs  such  as  MyoD,  Myogenin  (MyoG)  and  MEF2 
towards  various  muscle-specific  genes  and  promote  myogenic  differentiation.  Intriguingly, 
during  the  middle  period,  enhancer  sites  associated  with  enriched  MyoD  or  MyoG  binding  did 
not  exhibit  many  overlapping  enriched  motifs.  MyoG  peaks  at  enhancers  showed  enrichment  for 
multiple  ETS-domain  proteins,  Teads,  IRFs,  MEF2a  and  Stats,  which  was  in  contrast  to  MyoD 
peaks  at  enhancers  that  showed  enrichment  for  RunX,  Pbx,  p300,  Meisl  and  NF-k(3,  which  is 
modulated  by  the  PI3k-Akt  signaling  pathway.  Amalgamating  these  results  suggests  the  onset  of 
myoblast  differentiation  and  myocyte  formation  and  a  unique  cooperativity  for  different 
myogenic  TFs  (Fig.  5d). 

Positive  and  Negative  Regulation  of  Myogenic  Differentiation  Occurs  on  Cis  and  Trans  Levels 

Skeletal  muscle  regeneration  is  tightly  regulated  through  a  variety  of  positive  and  negative 
feedback  regulatory  pathways  and  an  excellent  example  is  through  the  IGF  pathway,  which 
serially  activates  PI3K,  Akt  and  mammalian  target  of  rapamycin  (mTOR).  Several  of  the 
effector  kinases  of  these  pathways  (protein  kinase  B  -Akt/PKB  and  extracellular-signal-regulated 
kinase  1/2  -  ERK1/2)  were  shown  in  Figure  3b  as  key  hubs  within  the  transcriptional  network 
orchestrating  the  early  phase  of  repair  and  regenerative  response.  One  component  of  the  IGF 
pathway,  the  PI3K  component,  is  essential  for  skeletal  muscle  regeneration;  and  the  PI3K 
enzymes  are  heterodimers  composed  of  a  catalytic  subunit  bound  to  a  regulatory  subunit.  Three 
Class  IA  PI3K  pi  10  catalytic  subunits  (a,  (3.  and  5)  are  expressed  in  skeletal  muscle34,  and  of 
these,  pi  10a  and  pllOp  have  been  shown  to  positively  influence  myoblast  proliferation  and 
differentiation35,36.  The  single  Class  IB  PI3K  pi  10  catalytic  subunit  y  is  expressed  abundantly  in 
immune  cells,  and  association  with  GPy  subunits  mediates  downstream  signal  transduction.  In 
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the  early  and  middle  phases  (24-72h)  following  injury  (the  period  associated  with  peak  myoblast 
proliferation),  we  observed  differential  enhancer  binding  and  increases  in  expression  of  PI3K 
catalytic  subunits  pi  10a  ( Pik3ca ),  p  I  1 0(3  ( Pik3cb ),  and  pi  105  ( Pik3cd)  for  the  injured  muscle  as 
compared  to  the  contralateral  control.  Figure  6  shows  increases  in  expression  and  active 
enhancers  of  the  PI3K  p85a  ( Pik3rl )  regulatory  subunit  along  with  increases  in  expression  and 
chromatin  remodeling  of  Class  IB  pllOy  ( Pik3cg )  and  its  regulatory  subunit  plOl  ( Pik3r5 ) 
during  the  healing  process. 

The  temporal  increases  in  PI3K  expression  (and  other  positive  myogenic  regulatory  pathways) 
between  24h-72h  post-injury  were  also  mirrored  by  chromatin  changes  and  increases  in 
expression  of  negative  feedback  regulators  such  as  MG53  (Trim7237,  XBP138,  USF1,  Id  1 ,  Id2, 
NFATcl39  and  several  HDACs),  which  associate  with  MyoD  and  MEF2  proteins  to  repress  their 
activity40.  The  negative  regulators  began  to  decrease  in  expression  starting  at  72h,  further 
reinforcing  the  observation  of  induction  of  myogenic  differentiation  at  this  time.  Collectively, 
these  observations  suggest  that  the  IGF,  PI3K  and  AKT  pathways  crosstalk  with  negative 
regulators  at  different  levels  to  direct  competition,  co-occupancy  or  occlusion  of  MyoD  and 
MEF2  binding  sites  (XBP1,  USF1),  binding  to  MyoD  and  MEF2  to  reduce  their  binding  affinity 
(Id  1 ,  Id2,  NFATcl,  HDACs)  as  well  as  disrupting  upstream  IGF  signaling  via  ubiquitin  ligases 
(MG53)  to  precisely  tune  myogenic  repair  and  regeneration  programs. 

Positive  and  negative  regulation  of  muscle  regeneration  has  also  recently  been  shown  to  occur 
post-transcriptionally,  where  muscle-specific  microRNAs  (miR-l-a/133a-2,  miR-l-b/133a-l, 
miR206/133b)  bind  to  the  3'  untranslated  region  (UTR)  of  mRNAs  that  promote  or  repress 
myogenic  gene  expression  programs  and  inhibit  translation41'42.  Small  RNA-Seq  was  perfonned 
at  multiple  time  points  from  each  stage  (early:  3h  &  lOh;  middle:  72h,  168h,  336h;  late:  504h, 
672h)  to  probe  trans -regulation  of  muscle  regeneration  dynamics  and  841  miRNAs  were 
detected  for  at  least  one  time  point  and  of  these,  143  miRNAs  showed  dynamic  behavior  (See 
Methods).  Figure  7a  shows  the  results  of  the  small  RNA-Seq  data  and  examination  of  muscle- 
specific  microRNAs  revealed  changes  in  expression  in  the  middle  and  late  periods 
(downregulation  of  miR-133a  and  miR-lb,  upregulation  of  miR-206,  miR-133b  and  miR-la), 
which  is  consistent  with  the  temporal  upregulation  of  their  regulatory  factors  MyoD  and  MyoG 
and  chromatin  profiles  exhibiting  active  enhancer  activity  (Fig.  7b)43.  Further  analysis  of  the 
datasets  reinforced  many  of  the  results  observed  from  the  RNA-Seq  and  ChIP-Seq  datasets, 
where  early  upregulated  microRNAs  were  associated  with  inflammation  and  immune  system 
programs  and  middle  and  late  stage  upregulated  microRNAs  were  associated  with  muscle  repair 
and  regeneration.  For  example,  in  the  early  period  the  inflammatory  miR-223  was  upregulated 
as  to  where  its  expression  subsided  in  the  middle  and  late  periods.  In  the  middle  period,  when  the 
nascent  stage  of  myogenic  differentiation  was  observed,  upregulation  of  miR126b  (associated 
with  angiogenic  and  chemokine  signaling44),  and  miR-429  (transcriptional  repressor  of  c-myc 
and  proliferation),  was  detected.  Additionally,  several  other  miRs  linked  to  skeletal  muscle 
regeneration46’46  exhibited  changes  in  expression  during  the  middle  period  such  as  the  miR-29 
family  (regulator  of  fibrosis  and  collagen  expression),  miR-2 1  (regulator  of  PI3k/Akt  signaling), 
miR-31  and  as  well  as  the  miR378  cluster  (regulator  of  mitochondrial  metabolism  and  inhibitor 
of  musculin).  The  delicate  balances  in  expression  of  these  /ra/?.s-regulatory  elements  through 
time  are  uniquely  further  titrated  by  competing  endogenous  RNAs  (ceRNAs)  such  as  linc- 
MD147,  which  sequesters  transcripts  such  as  miR133  to  facilitate  myogenic  differentiation 
transcriptional  programs.  Figure  7c  demonstrates  the  temporal  expression  profile  of  linc-MDl, 
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which  was  dramatically  upregulated  beginning  at  72h  and  peaked  at  168h,  consistent  with  the 
peak  times  observed  for  myoblast  proliferation  and  onset  of  differentiation. 

Discussion 

Muscle  recovery  after  traumatic  injuries  such  as  an  occupational  crush  or  blast  suffered  during 
military  combat  canonically  induces  a  predisposition  for  additional  injuries  and  chronic  pain. 
Many  of  the  various  muscle  repair  and  regeneration  sub-processes  triggered  after  trauma  have 
been  studied  and  shown  to  be  orchestrated  by  gene  expression  networks  controlled  by  lineage- 
specific  TFs  such  as  MyoD,  MyoG,  and  MEF2.  However,  many  questions  at  the  molecular  level 
remain  and  the  in-vivo  epigenomic  landscape  after  trauma  has  not  been  fully  characterized  and  as 
such,  our  understanding  of  how  these  different  factors  translate  to  influence  the  chromatin 
architecture  that  regulates  expression  at  lineage-specific  genes  has  been  limited.  Herein,  we  use 
integrative  genomic  mapping  technologies  to  profile  the  in-vivo  chromatin  state  of  various  cis- 
regulatory  elements  and  gene  expression  and  found  successive  waves  of  chromatin  remodeling 
orchestrated  by  combinations  of  lineage- specific  TFs. 

The  ChIP-Seq  maps  revealed  ds-regulatory  elements  such  as  promoters  were  largely  invariant 
when  compared  to  regions  demarcated  by  enhancer  marks  (H3K4mel,  H3K27ac),  which 
demonstrated  dynamic  behavior.  The  motifs  enriched  within  H3K27ac  regions  immediately 
after  the  injury  until  48-72h  later  were  associated  with  canonical  mediators  of  early  stress, 
immunity  and  growth  factor  responses  such  as  AP-1  (Fos/Jun),  SRF,  NF-k[1  EGR,  and  STATs. 
A  significant  fraction  of  the  TF  binding  sites  correlated  with  transcriptional  dynamics  and  was 
most  likely  attributable  to  infiltrating  monocytes  that  secrete  and  respond  to  cytokines  and 
chemokines.  These  early  immune-related  changes  were  paralleled  by  pathways  associated  with 
activation  and  proliferation  of  satellite  cells,  such  as  Notch,  Hippo  and  TGF-|3  signaling.  Three 
days  after  the  injury,  the  epigenetic  landscape  changed  and  could  be  attributed  to  proliferating 
satellite  cells  that  began  the  process  of  myogenic  differentiation.  These  changes  were 
orchestrated  through  concerted  binding  of  TFs  such  as  MEF2,  MyoD,  MyoG,  Teads,  RunX,  Pbx 
and  PPARy  and  chromatin  regulators48  such  as  p300,  CBP  and  SWI/SNF.  In  contrast  to  the  early 
period,  which  was  dominated  by  interactions  with  invading  immune  cells,  the  middle  and  late 
periods  were  enriched  for  TF-binding  sites  that  were  associated  with  structural  genes,  cellular 
migration  and  niche  remodeling.  The  two  different  types  of  interactions  offer  a  unique  insight 
into  how  the  plastic  chromatin  state  of  regenerating  muscle  changes  through  discrete  stages  and 
how  different  myogenic  TFs  require  cooperativity  to  execute  various  stage-specific  transitions. 

Intriguingly,  the  different  stages  of  healing  and  changes  in  chromatin  landscape  were 
orchestrated  only  when  binding  of  small  sets  of  lineage-specific  TFs49  that  co-localized  and 
cooperated  (or  antagonized)  with  other  TFs50  occurred.  An  example  of  this  titration  of  stage- 
specific  transitions  was  observed  between  24h-72h  post-injury,  whereby  temporal  increases  in 
members  of  the  IGF  and  PI3K  pathway  acted  to  induce  MyoD  binding  and  myogenic 
transcription  but  were  mediated  by  MG53,  which  augments  insulin  signaling,  and  XBP1,  USF1 
NFATcl  and  HDACs,  which  compete  for  E-box  binding  with  MyoD  and  MEF2  proteins.  This 
binding  pattern  uniquely  suggests  that  one  or  several  of  these  TFs  may  act  as  scaffolds  to 
cooperatively  recruit  various  other  TFs  and  chromatin  remodeling  complexes  to  enhancers  and 
tune  expression.  Cooperative  or  antagonistic  binding  of  TFs  during  different  stages  of  healing 
may  thus  provide  an  effective  means  of  reinforcing  or  tuning  expression  patterns51  and  facilitate 
a  metabolically  efficient  mechanism  to  quickly  respond  after  injury.  Since  muscle  repair  and 
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regeneration  utilizes  many  feed-forward  and  feedback  loops,  co-occupancy  or  competition  for 
binding  at  different  cis  elements  by  diverse  TFs  may  also  facilitate  a  precise  way  to  prevent 
extrinsic  signal  propagation  from  adjacent  tissues  that  are  also  regenerating  or  responding  to  the 
injury  and  rapidly  adjust  to  various  changes  in  metabolic  flow. 

The  temporal  events  of  inflammation,  regeneration  and  muscle  differentiation  are  also  fully 
recapitulated  at  miRNA  levels.  Disruption  of  the  muscle  integrity  stimulated  by  acute  trauma 
produced  miRNA  dynamics  very  similar  to  those  observed  in  other  muscle  myophathies, 
suggesting  the  modulation  of  common  molecular  pathways46'52. 

In  summary,  this  study  represents  one  of  the  first  studies  to  monitor  the  in-vivo  dynamics  of 
muscle  regeneration  at  genome- wide  levels.  The  simultaneous  measurement  of  epigenetic  marks 
and  RNA  sequencing  including  coding  and  noncoding  RNA  recapitulated  some  of  the  newest 
players  in  muscle  biology  and  allowed  unbiased  views  of  the  actions  of  various  TFs  on  their 
targets.  We  envision  that  a  wider  combinatorial  interrogation  of  such  a  dataset  can  represent  a 
valuable  resource  to  extend  the  networks  acting  in  such  a  complex  micro-environment  like  the 
cis-regulatory  modules  engaged  by  TF,  miRNAs  and  lincRNAs. 

Materials  and  Methods 

Animals  &  Traumatic  Injury  Model. 

Male  C57BL/6J  mice  (10  weeks  of  age,  24-27  grams)  were  obtained  from  The  Jackson 
Laboratory  (Bar  Harbor,  ME).  Mice  were  housed  one  per  cage  (shoebox  cage,  7"  x  1 1"  x  5"h)  in 
the  USARIEM  animal  facility  at  a  constant  Ta  =  24±l°C,  50%  relative  humidity,  with  a  12h / 
12h  (0600-1800  h)  light/dark  cycle.  Standard  laboratory  rodent  chow  and  water  were  provided 
ad  libitum.  Cages  were  supplied  with  Alpha-dri  and  cob  blend  bedding  for  nesting  and 
enrichment  and  plastic  houses  for  warmth  and  comfort.  Food  intake  and  body  mass  were 
recorded  daily.  Mice  were  cared  for  in  accordance  with  the  Guide  for  the  Care  and  Use  of 
Laboratory  Animals  in  a  facility  accredited  by  the  Association  for  the  Assessment  and 
Accreditation  of  Laboratory  Animal  Care  (AAALAC). 

Prior  to  administration  of  the  freeze  injury,  mice  were  anesthetized  with  a  combination  of 
fentanyl  (0.33  mg/kg),  droperidol  (16.7mg/kg),  and  diazepam  (5  mg/kg).  The  TA  muscle  was 
exposed  via  a  1  cm  long  incision  in  the  aseptically  prepared  skin  overlying  the  TA  muscle. 
Freeze  injury  was  performed  in  the  left,  hind  limb.  The  non-injured  contralateral  leg  served  as 
one  control.  Freeze  injury  was  induced  by  applying  a  6  mm  diameter  steel  probe  (cooled  to  the 
temperature  of  dry  ice,  -70C)  to  the  belly  of  the  TA  muscle  (directly  below  incision  site)  for  10 
seconds.  Following  injury,  the  skin  incision  was  closed  using  6-0  plain  gut  absorbable  suture 
(Ethicon,  Piscataway,  NJ).  The  analgesic,  Buprenorphine  (0.1  mg/kg  SQ)  was  administered 
using  a  25-27  gauge  needle  prior  to  recovery  from  anesthesia. 

Mice  were  euthanized  at  each  time-point  post-injury  (3,  10,  24,  48,  72,  168,  336,  504,  672  h)  via 
CO2  inhalation  (2  liters/min),  thoracotomy  and  exsanguination.  TA  muscles  were  removed  from 
the  injured  and  contralateral  limb;  weighed,  and  a  portion  of  the  tissue  was  crosslinked  in  1% 
formaldehyde  for  15  min  at  room  temperature,  with  constant  agitation.  The  reaction  was  then 
quenched  with  125mM  Glycine  for  5  min  at  room  temperature  with  constant  agitation.  The  fixed 
tissue  was  then  rinsed  3x  with  cold  PBS  (4°C),  spun  down,  snap  frozen  in  liquid  nitrogen  and 
stored  for  future  use. 

Chromatin  Isolation  and  Sequencing  Library  Preparation. 
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Each  frozen  tissue  was  thawed  for  at  least  30  min,  homogenized  (Tissue  Ruptor,  Qiagen)  and  the 
cells  were  lysed  for  at  least  10  min.  Next,  nuclei  were  isolated,  re-suspended  and  lysed  for  at 
least  10  min.  The  chromatin  was  then  sheared  using  a  Branson  sonificr  into  fragments  of  size 
range  of  150-700  base-pairs.  The  sheared  chromatin  was  then  incubated  overnight  at  4°C  with 
validated  antibodies  (Modified  Histone  Peptide  Array  -  Active  Motif)  H3K4me3  -  Millipore 
Catalog  #  07-743,  H3K27ac  -  Active  Motif  Catalog  #  39133,  H3K4mel  -  Abeam  Catalog  # 
ab8895)  using  constant  agitation  and  co-immunoprecipitated  using  a  mixture  of  Protein  A  & 
Protein  G  Dynabeads  (Life  Technologies)  for  2  hrs  at  4°C  with  constant  agitation.  The 
immunoprecipitated  chromatin  was  then  washed,  cross-links  reversed  and  treated  with  Proteinase 
K  and  DNA  was  purified.  Approximately  lOng  of  isolated  DNA  was  end-repaired  (End-It  DNA 
End-Repair  Kit,  Epicenter),  extended  and  A-tailed  (New  England  Biolabs),  and  ligated  to 
sequencing  adaptor  oligos  (Illumina).  The  adaptor-modified  library  was  then  amplified  by  PFU 
Ultra  II  Hotstart  Master  Mix  (Agilent)  and  size-selected  to  a  range  of  300-600  base-pairs  prior  to 
sequencing.  Libraries  were  pooled  and  sequenced  using  an  Illumina  Genome  Analyzer  IIx  using 
44  base-pair  single  end  reads  to  achieve  approximately  20  million  aligned  reads  per  sample. 

Quantitative  PCR  and  Small-RNA  Sequencing  Library  Preparation. 

A  portion  of  the  extracted  tissue  was  snap  frozen  in  Trizol  and  total  RNA  was  isolated  from 
the  tissue  after  homogenization  using  the  miRNeasy  Mini  Kit  (Qiagen)  as  per  the 
manufacturer’s  instructions.  RNA  concentration  and  integrity  were  measured  with  a 
Nanodrop  spectrophotometer  (Nanodrop  2000c)  and  Bioanalyzer  (Agilent  2100).  If  a 
sample  did  not  pass  quality  metrics  for  further  processing  (RIN>7),  the  samples  were  omitted 
from  further  processing.  mRNA  was  isolated  from  the  extracted  tissues  using  oligo(dT) 
Dynabeads  (Life  Technologies)  and  reverse  transcribed  to  cDNA  using  the  high  Capacity 
cDNA  Reverse  Transcription  kit  (Catalog  #  4374966;  Life  Technologies,  Carlsbad,  CA). 
TaqMan  MicroRNA  Reverse  Transcription  Kit  (Catalog  #  4366596)  was  used  to  generate 
template  for  microRNA  quantitative  PCR.  Quantitative  PCR  was  performed  using  conditions 
described  previously34  using  the  comparative  dCt  method.  Primers/probes  used  in  real-time 
PCR  are  shown  in  Supp.  Table  1. 

At  least  500  ng  of  isolated  total  RNA  was  used  to  produce  small-RNA  sequencing  libraries. 
Libraries  were  prepared  according  to  the  specifications  of  the  Truseq  smallRNA  kit  (Illumina) 
and  sequenced  on  a  MiSeq  single  ended  35bp  run  (Illumina).  A  custom  program  was  used  to 
match  the  FASTQ  sequences  to  the  microRNA  database  (miRBase,  www.mirbase.org)  for  the 
mouse  species  (mm9).  The  program  counts  exact  matches  and  matches  with  17+  base  identities 
with  two  or  fewer  mismatches  as  reads.  The  program  then  tallies  the  reads  for  each  miRNA 
detected  and  DESeq  was  used  to  normalize  the  miRNA  counts.  Heatmap.2  (R  language;  gplots 
module)  was  used  to  create  the  heatmap.  841  miRNAs  were  detected  with  RPKM>1  for  at  least 
one  time  point  and  143  miRNAs  showed  dynamic  behavior  (RPKM>100  for  at  least  one  time 
point  and  fold  change  >3  compared  to  uninjured  controls). 

ChIP-Seq  Data  Processing. 

Individual  replicates  were  aligned  to  the  mm9  mouse  genome  using  Bowtie253.  Replicates  for  a 
single  mark  at  a  single  timepoint  were  merged  with  the  samtools54  cat  command.  MACS255  was 
used  to  call  peaks  for  individual  replicates  (versus  the  merged  controls  for  a  given  timepoint)  as 
well  as  for  the  merged  datasets  (versus  the  merged  controls  for  a  given  timepoint).  The  fold 
change  parameter  in  MACS2  was  set  to  the  [2,500]  range,  and  the  -broadpeaks  option  was 
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utilized  to  call  peaks  for  the  K4mel  datasets.  The  -to-large  parameter  was  used  to  scale  the 
smaller  dataset  to  the  larger  dataset  and  the  extension  size  was  set  to  200  for  the  K4mel  datasets. 

The  standard  workflow  for  the  IDR  framework56  was  used  to  identify  irreproducible  peaks  from 
individual  replicates:  any  peaks  that  were  not  present  in  the  majority  of  replicates  were  excluded 
from  the  set  of  peaks  identified  by  MACS2  for  the  merged  datasets.  The  peaks  that  remained 
after  the  IDR  filtering  were  analyzed  with  the  GREAT  toolkit57  for  the  mm9  genome  to  identify 
genes  and  pathways  associated  with  the  enriched  sites.  The  gene  regulatory  domain  definition 
was  used  to  set  the  association  rules.  For  the  K4me3  marks,  the  GREAT  association  rules  were 
set  to  analyze  proximal  regions  1.5  kbp  upstream,  1.5  kbp  downstream,  and  distal  1.5  kbp.  For 
the  remaining  two  marks,  the  association  rules  were  set  to  analyze  proximal  regions  5.0  kbp 
upstream,  downstream,  and  distally.  The  MSigDB  pathways  identified  by  GREAT  with  FDR 
<0.05,  using  the  Bonferroni  correction,  were  selected.  Pathways  with  significant  FDR  values  at  8 
timepoints  for  K27ac,  7  timepoints  for  K4me3,  and  5  timepoints  for  K4mel  were  identified,  and 
significance  levels  were  plotted  over  time. 

In  a  parallel  analysis,  the  post-IDR  peaks  from  the  merged  replicates  were  compared  to  the 
genes. gtf  file  for  mm9.  Peaks  were  mapped  to  promoter  regions  (within  1.5  kbp  of  the  TSS),  the 
first  exon,  subsequent  exons,  and  introns,  as  determined  from  the  exons. gtf  file.  Peaks  were 
additionally  mapped  to  the  near-gene  region,  defined  as  within  5  kbp  upstream  of  the  TSS  or  5 
kbp  downstream  from  the  final  exon.  Remaining  peaks  were  assigned  to  the  intergenic  region. 

The  HOMER  software  suite58  was  used  to  discover  known  and  de  novo  motifs  for  the  K27ac, 
K4me3,  and  K4mel  datasets.  The  findMotifsGenome.pl  tool  within  the  HOMER  toolkit  was 
used  with  the  mm9  reference  genome,  and  a  threshold  of  le-50  was  used  to  determine  result 
significance  for  de  novo  motif  discovery. 

The  final  set  of  peaks  from  the  merged  datasets  were  compared  against  the  MyoD  ChIP-Seq 
datasets  from  Cao  et  al14  (GSM857390,  GSM857391)  as  well  as  a  MyoG  ChIP-Seq  dataset  from 
Wold  et  al15  (ENCFF001XVD).  The  bedtools  intersection  tool  was  used  to  identify  overlapping 
peaks  across  these  datasets  and  the  datasets  generated  as  part  of  this  study.  Upregulated  promoter 
regions  were  confirmed  by  analyzing  the  overlap  between  the  Myod/Myog  datasets  and  the 
K4me3  dataset.  The  K4me3  dataset  was  then  compared  against  the  K27ac  and  K4mel  datasets  to 
identify  and  remove  common  peaks,  leaving  only  K27ac  and  K4mel  peaks  that  map  to  enhancer 
regions.  These  remaining  peaks  were  subsequently  also  overlapped  with  the  MyoD/MyoG 
datasets  to  confirm  upregulated  enhancer  regions. 

The  KEGG  and  GO59  database  was  queried  to  identify  pathways  that  were  over-expressed  in  the 
injured  samples.  Additionally,  Gene  Set  Enrichment  Analysis60  using  the  MsigDB  v.  4.0  was 
performed  to  identify  differentially  enriched  gene  sets.  Using  the  Cytoscape  data  visualization 
software61,  upregulated  immune  response  pathways  and  gene  sets  (cytokine-cytokine  receptor 
inteaction,  chemokine  signaling  pathways,  leukocyte  transendothelial  migration,  genes  involved 
in  platelet  activation  and  signal  aggregation,  Jak-Stat  signaling  network,  AP-1  transcription 
factor  network,  genes  involved  in  signaling  by  interleukins,  118-and  Cxcr2 -mediated  signaling 
events,  HIF-1 -alpha  transcription  network,  TNF  receptor  signaling  pathway,  apoptosis,  NF-kB 
signaling  pathway)  and  GPCR  signaling  pathways  were  merged  to  construct  interaction  networks 
for  upregulated  genes  in  the  datasets.  All  member  genes  in  these  pathways  that  were  upregulated 
with  a  fold  change  of  two  or  greater  for  at  least  one  timepoint  were  included  in  the  interaction 
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network.  Gene  nodes  were  ranked  by  out-degree,  the  number  of  other  genes  in  the  network  that 
the  given  gene  interacts  with62. 

Data  Visualization 

The  IGV  browser  (v.2.3.36)63  was  used  to  visualize  pileup  data  and  to  generate  the  tracks  in 
Figures  1-5.  All  chromatin  maps  were  converted  and  visualized  as  bigWig  files  normalized  to 
10M  reads.  The  violin  plots  in  Supplemental  Figure  1  were  generated  with  MATLAB’s  violin 
function.  The  ggplots  heatmap2  function  in  R  was  used  to  generate  all  heatmaps. 
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Figure  1.  Experimental  overview  for  profiling  molecular  mechanisms  governing  in-vivo  tibialis  anterior  (TA) 
muscle  regeneration  after  severe  trauma,  a)  Schematic  diagram  of  injury  model  and  process  flow  for  chromatin 
and  mRNA  extraction  and  analysis  of  genomic  maps.  In  I,  the  TA  muscle  of  a  C57BL/6J  mouse  is  administered  a 
freeze  injury  and  the  contralateral  TA  muscle  serves  as  the  control.  The  injured  and  uninjured  tissues  are  extracted  at 
different  times  after  the  injury  was  administered  (3hrs-672  hrs).  During  II,  the  tissue  is  digested  and  mRNA  and 
chromatin  are  extracted  and  enriched.  In  III,  the  enriched  nucleic  acids  are  then  sequenced,  aligned  to  the  mm9 
genome  and  analyzed.  A  representative  example  of  the  Mef2a  gene  at  3  hours  post-injury  is  shown  where  the 
promoter  (labeled  P  in  gray)  and  enhancer  regions  (labeled  El  and  E2  in  purple)  are  depicted,  b)  Bar  graphs  of  gene 
expression  values  of  six  different  genes  corresponding  to  different  stages  of  the  muscle  regeneration  process  through 
time  from  left  to  right. 
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Figure  2.  Distribution  of  enriched  sites  across  the  genome  during  in-vivo  muscle  regeneration,  a)  Total  and 
differential  number  of  enriched  sites  for  H3K27ac  and  H3K4me3.  b)  Distribution  of  peaks  across  various  genomic 
elements  (promoter,  intergenic,  intron,  near-gene:  ±  2kbp  from  transcriptional  start  site,  exon),  c)  Shared  and  unique 
numbers  of  enriched  chromatin  sites  corresponding  to  different  genomic  elements  in  b).  d)  Representative  motifs  of 
overrepresented  transcription  factors  (p  <  10"10)  enriched  from  H3K27ac  and  H3K4mel -marked  regions  for  the 
different  stages  (early,  middle  and  late). 
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Figure  3.  Chromatin  landscapes  are  immediately  modified  after  severe  muscle  trauma  and  reflect  immune 

cell  infiltration,  a)  Heatmap  of  overrepresented  immune  system  pathways  derived  from  enriched  H3K27ac  peaks. 

b)  Interaction  network  diagram  of  upregulated  gene  sets  derived  from  over-represented  pathways  in  the  early  period. 

c)  Normalized  H3K4mel,  H3K27ac,  and  H3K4me3  ChIP-Seq  profiles  of  peaks  found  around  the  c-Fos  locus. 
Enriched  enhancer  regions  are  highlighted  in  gray,  d)  Bar  graphs  of  individual  gene  expression  values  of  four 
regulator  genes  through  time  from  left  to  right  (injured  samples  are  colored  and  uninjured  samples  are  uncolored). 
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Figure  4.  Injured  muscle  microenvironment  provides  cues  for  chromatin  remodeling  at  genes  associated  with 
satellite  cell  migration  and  proliferation,  a)  Line  graphs  of  individual  gene  expression  values  through  time  from 
left  to  right  of  SMAD-activated  genes  (TGF-|3  1  and  Tead4)  that  promote  TGF-|3  signaling  and  a  transition  towards 
Hippo  and  Wnt  signaling.  Serum  response  factor  (SRF)  targets  interleukin-4  (IL-4ra)  and  interleukin- 10  (IL-lOra) 
also  show  upregulation  at  this  time,  which  induce  macrophage  M2  polarization  and  modulate  MyoD  expression,  b) 
Heatmap  of  overrepresented  pathways  corresponding  to  G-protein  coupled  receptors  and  Rho  family  of  small 
GTPases,  which  were  derived  from  enriched  H3K27ac  peaks,  c)  Normalized  ChIP-Seq  tracks  of  H3K4mel, 
H3K27ac,  and  H3K4me3  profiles  around  the  myoferlin  (MyoF)  gene,  a  plasma  membrane  protein.  Enriched 
enhancer  regions  are  highlighted  in  gray  and  corresponding  enriched  TF  motif  is  labeled,  d)  Bar  graphs  of  individual 
gene  expression  values  of  MyoF  and  associated  transcription  factor  (TGIF1)  through  time  from  left  to  right  (injured 
samples  are  colored  and  uninjured  samples  are  uncolored). 
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Figure  5.  Coordinated  activation  of  myogenic  differentiation  and  healing,  a)  Analysis  of  total  and  differential 
activity  at  enhancers  for  MyoD  and  MyoG  derived  from  H3K27ac+  &  H3K4me3"  sites,  b)  Expression  heatmap  of 
IGF  signaling  genes  through  time  from  left  to  right,  c)  Enriched  GO  terms  derived  from  differential  chromatin  peaks 
at  enhancer  elements  show  associations  with  myogenic  differentiation.  The  corresponding  false  discovery  rates  for 
different  times  after  injury  are  plotted  and  the  size  of  each  circle  corresponds  to  the  number  of  significant  genes 
associated  with  the  category  (insulin  growth  factor:  IGF,  fibroblast  growth  factor:  FGF,  p38  mitogen-activated 
protein  kinase:  p38-MAPK,  phosphatidylinositol  3-kinase:  PI3K).  d)  Proposed  integrated  model  of  detected 
transcription  factors  (TFs)  associated  with  enhancer  elements  at  different  stages  after  injury. 
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Figure  6.  PI3K/Akt  activation  promotes  transition  from  proliferation  towards  myogenic  differentiation,  a) 

Bar  graphs  of  individual  gene  expression  values  for  various  components  of  the  PI3K  and  AKT  pathways  through 
time  from  left  to  right,  b)  Normalized  ChIP-Seq  tracks  of  H3K27ac  profiles  showing  differential  enhancer  activity  of 
the  PI3Krl,  PI3Kr5  and  PI3Kr6  loci.  Enriched  enhancer  regions  are  highlighted  in  gray  and  corresponding  enriched 
TF  motifs  are  labeled  underneath. 
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Figure  7.  Dynamics  of  non-coding  RNAs  after  severe  muscle  trauma,  a)  Expression  heatmap  for  small  RNAs  are 
plotted  through  time  from  left  to  right.  Muscle-specific  miRNAs  and  miRNAs  with  dynamic  expression  exhibited 
temporal  events  of  inflammation,  repair  and  regeneration.  Early  upregulated  microRNAs  were  associated  with 
inflammation  and  immune  system  programs  (miR-223)  and  middle  and  late  upregulated  microRNA  genes  were 
correlated  with  muscle  repair  and  regeneration  (miR-126b,  miR-29,  miR-21,  miR-31,  miR-378).  b)  Chromatin 
profiles  of  several  microRNA  genes  for  two  histone  modifications  (H3K27ac-blue  and  H3K4mel -orange) 
associated  with  enhancer  elements.  The  maps  are  ordered  through  time  and  MyoD  enrichments  from  Cao  et  al14  are 
labeled  underneath  the  track  with  a  pink  box.  Each  scale  bar  is  10  kbp.  c)  Expression  profile  of  the  linc-MDl  gene, 
which  sequesters  microRNAs  that  inhibit  transcriptional  programs  such  as  myogenic  differentiation.  The  expression 
of  linc-MDl  increased  at  72h  and  peaked  at  168h,  which  was  consistent  with  the  peak  times  observed  for  myoblast 
proliferation  and  onset  of  differentiation. 
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